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Abstract. We present a set of optimum ground states for a large class of 
spin- 1 chains. Such global ground states are simultaneously ground states of 
the local Hamiltonian, i.e. the nearest neighbour interaction in the present 
case. They are constructed in the form of a matrix product. We find three 
types of phases, namely a weak antiferromagnet, a weak ferromagnet, and a 
dimerized antiferromagnet. The main physical properties of these phases are 
calculated exactly by using a transfer matrix technique, in particular mag- 
netization and two spin correlations. Depending on the model parameters, 
they show a surprisingly rich structure. 
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1 Introduction 



Quantum spin systems have been investigated extensively in recent years, 
in particular in one dimension. Experimentally many quasi-one-dimensional 
materials have been discovered for quantum spin 5 = i,l,|,2,... and in- 
teresting physical properties have been observed [1-9]. Most experiments 
indicate that such systems have restricted symmetries and that anisotropy 
effects are quite important. 

Theoretical interest has been mostly initiated by Haldane [10,11] who 
emphasized essential differences between integral and half-integral spin sys- 
tems which could lead to qualitatively different physical properties. For 
integer spin he showed the existence of a new antiferromagnetic phase, i.e. 
the Haldane phase, characterized by a unique ground state, a gap to the ex- 
citations, and exponential decay of correlations. In subsequent work exact 
ground states have been constructed for such integer spin models, mostly for 
spin 1, and mathematical properties of such states have been investigated 
in much detail (see for instace [12-14]). 

In recent work [15-18] we have considered a most general class of aniso- 
tropic spin-1 chains. Exact ground states were found in a large parameter 
subspace of this model in the form of a matrix product and have been called 
matrix product ground states. In the present paper (also [|l^) we shall apply 
the same idea, namely the construction of exact ground states via matrix 
products, to the case of spin S = ^. It turns out that in this case three 
different types of homogeneous states can be constructed, namely a weak 
antiferromagnet, a weak ferromagnet, and a dimerized antiferromagnet, all 
of them characterized by a broken symmetry, and thus by long-range order. 
Of course we find no Haldane antiferromagnet. 

The paper is organized as follows. In section ^ we introduce the notion 
optimum ground state for a quite general lattice model whose global Hamil- 
tonian is the sum of suitably defined local Hamiltonians, for instance nearest 
neighbour interactions. It has the property that this global ground state is 
simultaneously ground state to all local Hamiltonians. 

In section ^ we consider and classify the Hamiltonian for most general 
spin- 1 lattice models with nearest neighbour interaction exhibiting some 
minimum set of symmetries. The parameter space turns out to be 12- 
dimensional. 

In section ^ it is shown that optimum ground states in the sense of section 
^ can also be constructed for spin-| chains in the form of matrix products 
provided certain conditions are satisfied. These will lead to restrictions on 
the interaction parameters of the general model in the following subsections. 

Section ^ contains the main results. In three subsections three different 
types of optimum ground states in the form of matrix products are pre- 
sented. As mentioned before they describe a weak antiferromagnet, a weak 
ferromagnet, and a dimerized antiferromagnet. In each case the single site 
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state matrices and the parameter restrictions are determined. Also the most 
important ground state properties are expHcitly calculated, such as magne- 
tizations, two-spin correlations, and corresponding correlation lengths. All 
truncated correlations show exponential decay. The main steps of calcu- 
lating such properties via transfer matrix technique are briefly derived in 
Appendix A. 

In section ^ we describe an alternative graphic representation of optimum 
ground states as a vertex state model. This is formulated for the spin-| chain, 
but can be generalized to other spin values S and to applications in higher 
dimensions [^0|. Section ^ contains a short summary of main results. 

2 Optimum ground states 

We consider a spin system on an arbitrary lattice with homogeneous nearest 
neighbour interaction. The Hamiltonian of such a system is of the form 



Angular brackets denote nearest neighbour lattice sites i and j, to 

which spin variables are attached. As the system is homogeneous, all local 
Hamiltonians hij are in fact equal, they only act on different spin pairs. By 
simply adding a trivial constant to the Hamiltonian - which never changes 
the physics of the system - one can always achieve that the lowest eigenvalue 
of h vanishes, i.e. eo = 0. In this case, the lowest eigenvalue of -ff - i. e. the 
global ground state energy Eq - is greater than or equal to zero, since the sum 
of non- negative operators is also non- negative. Now one has to distinguish 
between two cases: 
A) So > 

This is the usual case. The ground state of the complete spin system 
involves not only local ground states of h but also (possibly all) excited 
states of the local Hamiltonian. As there is no local condition for the global 
ground state, one is dealing with a true many body problem. Therefore the 
construction of the global ground state of such a system is usually extremely 
difficult or even impossible. 



In this very special case, there exists a local condition for finding the 
ground state. Let |^o) be the ground state of H. Then we have 



Since all hij are non-negative operators, every single (^o|/ii,j|^o) must van- 
ish. If A is the largest eigenvalue of hij, we have the estimate 




(1) 



B) = 




(2) 




(3) 
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This of course requires 



/ii,j|^o) = 



for all 



(4) 



This means: The global ground state consists only of ground states of the 
local Hamiltonian, no excited local states are involved. A ground state 
of this type is called optimum ground state. To obtain optimum ground 
states for a given system we have to perform two steps. First, the ground 
states of the local Hamiltonian must be determined, which is a local problem 
with few degrees of freedom. Second, one has to check, whether these local 
ground states can be stuck together to form a global ground state. The 
most important feature of optimum ground states is the existence of the 
local condition (^) . In section ^ we show how optimum ground states can 
easily be constructed by using matrix products. 



The physically most interesting systems are those exhibiting a certain degree 
of symmetry, as they are most often realized in experimental situations. 
We therefore restrict ourselves to Hamiltonians with only nearest neighbour 
interaction (for simplicity) and the following set of symmetries: 

• Homogeneity: All local Hamiltonians hij are equal, they only act on 
different spin pairs. 

• Parity invariance: /ij.j commutes with the parity operator Pij, which 
interchanges the two neighbour spins i and j. 

• Rotational invariance in the xy-plane of spin space: hij commutes 
with the pair magnetization operator Sf + Sj. 

• Spin flip invariance: — > —S^ leaves hij unchanged. 

In the following, we focus on spin-| systems. In this case, the operators 
hi J, Pi J and Sf + 5| act on a 16-dimensional vector space (the space of 
a spin pair). These operators commute with each other, so they can be 
diagonalized simultaneously. It turns out that Pij and Sf + 5| divide the 
vector space of a spin pair in small subspaces of dimensions less or equal to 2. 
Therefore, the eigenstates of the local Hamiltonian are restricted to not more 
than 2-dimensional invariant subspaces. This enables us to give a simple 
parametrization of all possible local Hamiltonians in terms of projectors 
onto the local eigenstates \vk)- 



3 Spin 



-| Hamiltonians 



16 




(5) 



k=l 
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The Afc are real parameters that can be chosen arbitrarily without breaking 
one of the above symmetries. 

Within a two-dimensional invariant subspace, the two eigenstates \vk) of 
h are not uniquely fixed by the operators Pij and Sf + Sj. Instead, they 
form an ortho-normal basis that may be rotated by a certain angle. Besides 
the A-parameters, one therefore has to introduce one additional parameter 
for each two-dimensional invariant subspace in order to take care of these 
possible rotations. 

A convenient notation for the canonical spin-| basis is 

sfm = ii3) sf\i) = 

(6) 

Sf\3) = -i|3) Si\l) = . 

The following table contains a classification of all eigenstate^ of h according 
to the corresponding eigenvalues n and p of the operators Sf + Sj and Pij, 
respectively: 



M = 


3, 


p = 


1 


: t;3 


= |33) 




M = 


-3, 


p = 


1 


: t;_3 


= |33) 




^J' = 


2, 


p = 


1 




= |31) + 


|13) 


M = 


2, 


p = 


-1 


: V2 


= |31)- 


|13) 


M = 


-2, 


p = 


1 


■ 


= |3T) + 


|13) 


M = 


-2, 


p = 


-1 


: v-_2 


= |3T)- 


|T3) 


= 


1, 


p = 


1 




= + 


f(|3T) + |T3)) 


M = 


1, 


p = 


1 


■ ^12 


= - 


-(|3T) + |T3)) 


fj- = 


1, 


p = 


-1 


■ 


= |3T)- 


|T3) 


^J' = 


-1, 


p = 


1 


■ ^lll 


= |TT> + 


f(|31) + |13)) 


fj- = 


-1, 


p = 


1 




= a\Tl) - 


-(|31) + |13)) 


M = 


-1, 


p = 


-1 


: ^^Zi 


= |31)- 


|13) 


M = 


0, 


p = 


1 


: ^01 


= (|11) + |11)) + 6(|33) + |33)) 


M = 


0, 


p = 


1 


: ^^02 


= Kill) 


+ |Tl))-(|33) + |33)) 


M = 


0, 


p = 


-1 


: ^01 


= (|1T)- 


-|Tl)) + c(|33)-|33)) 


M = 


0, 


p = 


-1 


: Vo2 


= c(|lT) 


-|T1))-(|33)-|33)) 



(7) 



The states are not normalized, because all normalization factors can be absorbed into 
the A-parameters. 
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States with p = 1 are symmetric under spin transposition, the others are 
antisymmetric. The parameters a, b, and c serve to handle the possible 
rotations in the two-dimensional invariant subspaces. 

Because of spin flip invariance, not all A-parameters are independent. 
A spin flip transforms an eigenstate with pair magnetization fi into a cor- 
responding eigenstate with pair magnetization — /i. These two eigenstates 
must have the same energy and therefore the same A-parameter. This leads 
to the following most general form for all local spin-| Hamiltonians which 
exhibit the required set of symmetries: 

h = X3{\v3){v3\ + \v^3){V-3\) + 

^t{\vt){vt\ + \vl2){V-2\) + 
>^2i\v2){^2 \ + l^'-2)("-2l ) + 

>^tii\vti){vti\ + \^^ii){v^ii\)+ 

A|2(kl2)Hl + l«^12)(«^12l) + 
\ + ) + 

^01 1 «01 ) (^'Ol I + ^02 1 ^02 ) (^^02 1 + 
-^01 1 "01 ) (^'Ol I + -^02 1 ){V02\ ■ 

The parameter space of all possible local Hamiltonians is 12-dimensional 
by fixing the local ground state energy at zero, i.e. cq = A^^jj^ = (as 
required in section |2|), as we are left with 9 A-parameters (> 0) and the 3 
superposition parameters a, 6, c. These 12 parameters include a trivial scale 
parameter, thus we have 11 non-trivial (interaction) parameters. 

We briefly mention the usual representation of the local interaction /ijj- 
in terms of the canonical spin operators 5^, 5^, 5^. Because of the required 
symmetries, only polynomials of 3 pair operators are allowed, 

U = SfSJ + SfS] 

V = SfS^ (9) 
W = (5f)2 + (5|)2 . 

Furthermore, for spin | at most second powers of spin operators are required, 
as higher powers are reduced by the spin algebra. Then it is easily seen that 
the most general interaction hij can be written as 

hij = const. + aiU + a2V+ 

asU^ + a^iUV + VU) + agF^ + ^^^^ 

arU^ + asUVU + agVUV + aiQV^+ ^ ' 

au{UW + WU) +ai2{VW + WV) , 
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in terms of 12 interaction parameters Of course, this representation 
(10) is completely equivalent to (^). We refrain from recalculating the a- 
parameters in terms of the A-parameters and a, h, c and vice versa. This is 
simple algebra and is not needed in the following. 



4 Matrix product representation 

In one dimension, i. e. for spin chains, optimum ground states can easily 
be constructed by using matrix products. This method was introduced in 
[15-17] and extended in jl^ to local Hamiltonians with next-nearest neigh- 
bour interactions. The Hamiltonian for a chain of length L with periodic 
boundary conditions is 

L 

H = Y,h,,,+i . (11) 
1=1 

Now consider a 2 x 2 matrix with single spin states as its elements 

m = . (12) 

V 1*3) \S4) J 

The product of two such matrices is defined as 

[in- g)ki = ^mkn® Qni , (13) 

n 

i. e. the tensorial product is applied to the elements. Obviously, the trace of 
an L-fold product of such matrices is a valid L-spin state and can serve as 
an ansatz for the groundstate of the system. In principle, all these matrices 
along the chain can be different, but in our case it is sufficient to use at most 
an alternating product of only two different matrices^ 

l^'o) = tr mi -52 • (14) 

To be an optimum ground state of the spin chain, |^'o) has to be annihilated 
by all local Hamiltonians This leads to the following theorem: 

If the local Hamiltonian annihilates all elements of the two ma- 
trix products m ■ g and g ■ m, then |^o) is an optimum ground 
state of the global Hamiltonian H. 



5 Ground states for spin-| chains 

In the following subsections we construct three different types of optimum 
ground states for spin-| chains using the matrix product technique. We also 
calculate the single spin magnetization and two point correlation functions. 

^To avoid possible misfits for closed chains we shall also restrict ourselves to even values 
of L. 
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5.1 Weak antiferromagnet 



An interesting set of ground states is obtained by using the alternating 
(along the chain) product of the two matrices 



m 




9 




(15) 



with the discrete parameter o" = ±1 and an arbitrary real parameter a. The 
elements of the two matrix products m ■ g and g ■ m are the pair states 



11) 


+ 


a|3T) 


11) 


+ 


aT3) 


11) 


+ 


a 31) 


TT) 


+ 


a 13) 


IT) 


+ 


cja2|33) 


11) 


+ 


aa^ 33) 


11) 


+ 


fjlTl) 



(16) 



If h annihilates these 7 pair states, then 
iV'o) = tr mi • 52 • m.3 • 54 • • • 



(17) 



is an optimum ground state of the global Hamiltonian H. This is achieved 
if the following constraints are imposed on the parameters defined in equa- 
tion (Isl): 



A 



11 - 
b = c 



a; 



a;^ 



1 — ^'01 
(TO? 

a arbitrary, a 



^02 



^01 







^12' ^02*^ > 
= ±1 . 



(18) 



The inequalities ensure, that all other pair states have a higher energy level. 
In this case, the local Hamiltonian h has 7 ground states and 9 excited 
states. The system still contains 6 continuous and 1 discrete parameter (a). 
Expectation values of observables can be calculated using the transfer 



matrix technique, which is explained in [15|, for any given chain length L. 
However, we shall only give the results for the thermodynamic limit L — > oo, 
since this is the most important case. The transfer matrix and its eigensys- 
tem can be found in Appendix |A.l| . The most interesting expectation value 
is the single spin magnetization. If {S^^g)^^ denotes the magnetization of a 
spin, represented by an m- or g-mathx, respectively, we have 



Ho 



a 



g/i>o 



1 



1 



ip 



_i_ 3 

75'2 



(19) 



which is calculated from the formulae of Appendix A. So the single spin mag- 
netization alternates from site to site, which indicates an antiferromagnet 
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with Neel ordering. The global magnetization vanishes. As the sublattices 
are only fully polarized if a ^ cxd, we call this ground state a weak antifer- 
romagnet. 

Unlike the single spin magnetization, the average of the square of the 
operator covers the full range of possible expectation values: 



1 



+ 



4a2 



1 9 
4' 4 



(20) 



4 (a2-l)2+4-{a2-3)^(a2-l)2+4 

Because of rotational invariance in the xy-plane, this also gives 



3/3 



2 V 2 



+ 1 



7 



2a2 



(21) 



4 (a2-l)2+4-(a2-3)V(a^-l)H4 

The next interesting observable is the nearest neighbour correlation 



z\2 _1 l-a^+V(a2-l)2+4 



(22) 



'm/^o 2 ■ (a2-l)2+4 ' 2+^{a^-ir+4 

which is plotted as a function of the parameter a in figure |l|. For large values 



i '-'j+l/'/'O 



Figure 1: Nearest neighbour correlation of the weak antiferromagnet 

of a, the nearest neighbour correlation approaches — |, which is the value 
for a strong antiferromagnet with fully polarized sublattices. 

The longitudinal correlation function consists of a constant and an ex- 
ponentially decaying term: 



{StS'r)^o = {S'Jlo + for odd 



(23) 
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where is the longitudinal correlation length. Its inverse is given by 



-1 



In 



1 + a2 + V(a^- 1)^ + 4 



1 + a2 - V(a2- 1)^ + 4 
and is plotted as a function of a in figure 



(24) 



Since the inverse correlation 




0.5 1 



Figure 2: Inverse longitudinal (thin) and transverse (thick) correlation 
length of the weak antiferromagnet 

length never vanishes, the system is non-critical. But at a = 1 it goes to 
infinity, i. e. the complete exponentially decaying term in the correlation 
function vanishes, only the constant term survives. This term indicates 
long-range antiferromagnetic order. 

Unlike the longitudinal one, the transverse correlation function (5f 
does not exhibit a constant term. It decays exponentially with inverse cor- 
relation length 



er' = In 



1 + a2 + V(«'-l)' + 4 



(25) 



which is also plotted as a function of a in figure It never vanishes nor 
diverges. 

Finally we note that in fact we have from (0) two ground states by 
exchanging sublattices, namely 



tr mi • 32 • "^3 • 54 • • • and 
tr 51 • m2 • 53 • "14 • • • . 



(26) 



These are the only ground states of the global Hamiltonian H for the closed 
chain. The proof can be carried out using a straightforward induction over 
the chain length L and is omitted here. This degeneracy indicates that the 
sublattice symmetry is broken, as is typical for a Neel antiferromagnet. 
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5.2 Weak ferromagnet 



Instead of using the alternating product of the matrices m and g defined in 
equation (^), it is also possible to construct global ground states from the 
homogeneous products 



{(Pq) = ti mi ■ 1712 ■ m^- ■■ rriL 
\iPq) = ti gi ■ 92 ■ 93- ■■ gL ■ 



(27) 



These two states are optimum ground states of H, if the local Hamiltonian 
h annihilates the elements of m • m and g ■ g, which are multiples of the 
following 7 pair states: 

31) + cj|13) 

31) + cj|13) 

11) + a|3T) 

11) + a|T3) (28) 

TT) + a|31) 

11) + a|13) 

11) + ct|11) . 

Therefore, if the parameters defined in equation (^) obey 



A? 



An - Aj 



Aqi 







6 = 

As; ^2'^^ A5*2, Aq]°", Aq2, ^q2 >0 

a, c arbitrary, 



(29) 



then \^Pq) and \^q) are optimum ground states of the global Hamiltonian H. 
The system still contains 8 continuous parameters (6 A-parameters and a,c), 
which are not fixed. The global ground state, however, only depends on one 
continuous (a) and one discrete parameter (a). 

As for the weak antiferromagnet, expectation values of local observables 
can be calculated using the transfer matrix technique. The eigenvalues and 



the corresponding eigenvectors can be found in Appendix A. 2, It is sufficient 
to calculate expectation values of \(Pq), since those of |(^q ) are obtained by 
replacing by —S^. The magnetization per spin is exactly 



(5f 



(30) 



independent of the parameter a, which indicates a state of ferromagnetic 
ordering. But unlike the strong spin-| ferromagnet, which has a magnetiza- 
tion of I per spin, the system is not fully polarized. We therefore call this 
state a weak ferromagnet. 

The fluctuations of Sf are o-dependent and given by 



{{St 



1 \a\ 
4 ^ iTR 



1 5 

4' 4 



(31) 
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Because of rotational invariance in the xy-plane, this also yields 



7 



4 2(l + |a|) 
The longitudinal two-spin correlation function is 



/ QZ QZ 



(i-H 



1 + \a\ 



for r > 2 



(32) 



(33) 



which consists of a constant term - indicating the long-range ferromagnetic 
order - and an exponentially decaying term with inverse correlation length 



= In 



1 + \a\ 



(34) 



It is plotted as a function of a in figure ^ Note the divergence at a = 1, 




Figure 3: Inverse longitudinal (thin) and transverse (thick) correlation 
length of the weak ferromagnet 



which also occured for the weak antiferromagnet and indicates complete 
absence of truncated correlations for r > 3. As special case we also get the 
nearest neighbour correlation 

1 

(5f = - - ^^^^ . (35) 

Its a-dependence is shown in figure |^. Although the global state has ferro- 
magnetic long range order, the nearest neighbour coupling is negative for 
\a\ > 1. This is due to the fact, that for large \a\ the two states with 
S"^ = I and = are dominating in the m-matrix ( p!5|) which leads to 
an alternation of these dominating states in the global state (p7|). 
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a 

Figure 4: Nearest neighbour correlation of the weak ferromagnet 

Unhke the longitudinal one, the transverse two-spin correlation function 
decays purely exponentially. The corresponding inverse correlation length 



One can prove by induction with respect to L = 2,3,4,..., that \(Pq) 
and Iv^o ) ^-re the only ground states of the global Hamiltonian H for the 
closed chain. The proof is completely analogous to the one for the weak 
antiferromagnet and is again omitted. 

5.3 Dimerized antiferromagnet 

Unlike the weak antiferromagnet and the weak ferromagnet discussed before, 
a last type of optimum ground states is not composed of square matrices. 
Instead, the alternating product of the 2x3/3x2 matrices, 





(37) 



can be used to construct the global ground state 




(38) 
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The elements of the products G ■ M and M ■ G are hnear combinations of 
the fohowing 9 pair states: 

1 31) + cj|13) 
|31) + cj|13) 
111) + a|3T) 
111) + a|T3) 

jll) + a|31) (39) 
|TT) + a|13)_ 
jll) + cja2|33) 
|Tl) + cja2|33) 
I IT) + cj|Tl) . 

These states are annihilated by the local Hamiltonian h if the parameters 
in equation dsl) satisfy 



A? 



'^n - AJ 



Aqi 



Ao2 







(40) 



b = c = aa 

As; A2 A]*2, Aq2'^ > 
a arbitrary, o" = ±1 

In this case, the local Hamiltonian has 9 ground states and 7 excited states. 
The system still contains 5 continuous and 1 discrete parameter (a). In the 
isotropic case, i. e. a = and a = —1, two spin correlation functions 

have already been calculated in |21, 22 1. 

The transfer matrix and its eigensystem can be found in Appendix |A.3| . 
Unlike the ground states presented in subsections 5.1 and |5.2| , this state is 
antiferromagnetic in the sense that (•S'^^otai) ~ ^ vanishing sublattice 

magnetization 

{Shho = {Shho = o . (41) 

As, however, the alternating G and M matrices are different, sublattice 
symmetry is broken. The global ground state is therefore a dimerized anti- 
ferromagnet. 

The fluctuations of are given by 

((•S'g')^)*o = (('S'm)^)*o = - - — — ^ — - ^ 



1 9 
4'4 



(42) 



4 (a2 + 1)2 + 

thus covering the complete range of possible values. Because of rotational 
invariance we also have 



4 (a2 + l)2 + ; 



(43) 



As G and M are not square matrices, there are two different nearest neigh- 
bour correlations in contrast to the situation in subsections 5.1 and If 
the left matrix is of G type, we have 

9 (a2 + 1)2 
'4'(a2 + l)2+8 ' 



('S'g"S'm)>I'o 



(44) 
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Figure 5: G ■ M (thin) and M ■ G (thick) nearest neighbour correlation of 
the dimerized antiferromagnet 

a, these two correlations do not coincide. The alternating binding structure 
characterizes the broken symmetry. This type of state is called dimerized. 

Both the longitudinal and the transverse correlation functions decay ex- 
ponentially 

(5f5^)vi>o = Ae-''/^' for odd values of r 
{SfS-r)s,, = Ate-^l^^ . 

The corresponding inverse correlation lengths 

1 _ V(a^+l)^+8 

- \ a'-l\ (47) 

= lniV(«' + l)' + 8 

are plotted as a function of a in figure ^. As both do not vanish for any 
finite value of a, there is no critical behaviour. At a = 1 the longitudinal 



correlation length again diverges as in subsections 5.1 and ^.21 . 

The list of local ground states ( |39D actually is the combination of all 
local ground states of the weak antiferromagnet (section |5.lD and the weak 



ferromagnet (section 5.2). Therefore it is clear that the corresponding global 



states are also optimum ground states of the global Hamiltonian with pa- 
rameters (^). However, it turns out that the degeneracy is even much 
higher. Any state of the form 

1^') = tr pi ■ P2 ■ P3 - ■ - PL (48) 
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a 



Figure 6: Inverse longitudinal (thin) and transverse (thick) correlation 
length of the dimerized antiferromagnet 



is an optimum ground state of the global Hamiltonian H, if every single 
matrix pi is one of the matrices m or g (|T^) from subsection 5^ . The reason 
for this is that the local Hamiltonian h annihilates the elements of all four 
possible products of m and g: 



h{m ■ m) = h[g ■ g) = h{m ■ g) = h{g ■ m) = 



(49) 



We therefore have a global degeneracy of at least no = 2^. Among these 
states are the following homogeneous special cases. The two states 



|<y3o) = tr mi • 7712 • 7713 ••• 

l<^o> = tr c/1 • c/2 • 53 • • • 



(50) 



represent the weak ferromagnet discussed in subsection 5.2. The two states 

HI) 



tr 7711 • 52 • • 54 • • • 
tr 51 • 7772 • 53 • "14 • • • 777L 



(51) 



represent the weak antiferromagnet discussed in subsection 5.1. Finally, we 
have the dimerized antiferromagnet (|38|) of this subsection. 



6 Vertex state model representation 

Apart from the foregoing matrix representation which is especially conve- 
nient for calculational purposes, it is also possible to give a graphic represen- 
tation of the constructed optimum ground states in terms of vertices. This 
vertex state model representation is especially useful to construct optimum 
ground states in higher dimensions pO|| . 
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For a specific matrix m we map each element on a vertex with two 
emanating bonds. A variable is attached to the left and the right bond 
corresponding to the left and right matrix index, respectively. The bond 
variables are represented by arrows pointing out of or into the vertex: 



mil 
mi2 
m2i 
m22 



< • < 
< • > 



(52) 



> • > 



The notation for a matrix product is 
(m • m)i2 ^ < ♦ < ♦ > 



->- 



+ < • > 



(53) 



i. e. a contraction denotes the sum over the bond variable. Since we are 
interested in matrices, whose elements are single spin states, we must asso- 
ciate a weighted spin state with each vertex rather than a number (vertex 
weight) in the usual case. The result will therefore be a "vertex state model" 
rather than a "vertex model". 

In order to describe global ground states composed of more than one type 
of matrices, it is also necessary to introduce different classes of vertices. This 
leads us to the following classification of vertex states corresponding to (111 



< • < 
< • > 



> • < 
> • > 



<Q> 
>Q< 
>Q> 



1) 

3> 
T) 

1) 

T) 
1) 
3) 
T) 



(m-class) 



(54) 



-class) 



As an example the following graph represents the ground state \ipo) of the 



weak antiferromagnet from subsection |5.l| for the closed chain with L 




(55) 
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We note that for a given class of vertices, the eigenvalue of the vertex 
state is given by 



= of outgoing arrows — # of ingoing arrows) 

+ 




(56) 

This also ensures conservation of the local pair magnetization Sf + Sfj^i- 



The 2x3- and 3 x 2-matrices from subsection 5.3 can also be handled 
easily using a new type of bond that carries a bond variable with values: 
'left arrow', 'zero', and 'right arrow'. 



7 Summary 

We have presented three different one-parameter sets of optimum ground 
states for the spin-| chain. The first one is a Neel ordered antiferromagnet 
which exhibits a parameter dependent sublattice magnetization. Pull polar- 
ization of the sublattices is only achieved in the limit a — > oo. We therefore 
call this state a weak antiferromagnet. It is an optimum ground state for a 
6-dimensional manifold of Hamiltonians with usual symmetries. The global 
degeneracy is 2. Truncated correlations decay exponentially. 

The second set of states is of ferromagnetic type. The global magneti- 
zation is which is one third of the maximum value. So this is a weak 
ferromagnet. The corresponding set of Hamiltonians is an 8-dimensional 
manifold. Again the global degeneracy is 2, corresponding to the broken 
—S^ symmetry. Also truncated correlations decay exponentially. 

Finally we have the dimerized antiferromagnet. Its properties are quite 
different from the first two types of ground states. There is no sublattice 
magnetization, but the broken sublattice symmetry here leads to an alternat- 
ing binding structure of nearest neighbour spins. The corresponding set of 
Hamiltonians is a 5-dimensional manifold, thus smaller than in the two other 
cases. Actually we have a vast number of ground states of the type (|4^), 
the degeneracy is at least 2^, corresponding to a specific zero-temperature 
entropy of In 2. 



Appendix A: Transfer matrices 



The transfer matrix method has been explained in detail in |15|. In the 
following, we only cite the main steps of this technique and summarize the 
results for the optimum ground states presented in section |5|. 

To calculate the expectation value {Ai) of a local operator Ai acting on 
spin states at site i within the state (jlj), we define the transfer matrix 



T(kuh),(k2,i2) = 9)kuk2\im- g)h,i2) ■ (A.l) 
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Redefining index pairs as 



(1,1) 
(2,2) 

(1,2) 

(2,1) 



1 

2 
3 
4 



(A.2) 



the tensor ( [A.l| ) results in a 4 x 4-matrix T. Additionally, for any operator 
A acting on (m • g) we define a tensor 



^(fci,h),(fc2,Z2)(^) = {{m- g)k^,k2 1^1 {m • 5)/i,/2 ) 



(A.3) 



which also reduces to a 4 x 4-matrix T{A). These two matrices can now be 
used to calculate the expectation value for a single spin operator Ai as 



(^ol^il^o) _ trace T{A)T^I^- 
(^'ol^o) ~ 



trace T^/^ 



(A.4) 



Usually, T is a symmetric matrix, so we can rewrite ( A.4 ) in terms of the 
eigenvalues Xn and eigenvectors Un of T: 



{A, 



L/2-1 
n 



L/2 



(A.5) 



Similar considerations also yield an expression for correlation functions. If 
Ai and Bj are operators acting on spins at sites i and j, respectively, we get 

\T{A)\Un2) X'n-,^ {Un,\T{B)\Un,) Xni' '\ , 

{AiBr) = — (^-6) 

l^n=l Xn 

as the final result for all correlation functions of two local operators. Because 
of translational invariance, we fix the first site at i = 1 and denote the second 
site j = r, r = 2, 3, . . . . 



A.l Weak antiferromagnet 

For the ground state I'i/'o) (P^, the transfer matrix ( [A.l| ) is given by 

\ 



1 
1/ 

The eigenvalues Xn and the corresponding eigenvectors m„ are 



( l + a 
l + a 


V 



4 l + a2 
2 2 





(A.7) 



Xl = {l + X + ^Jl + X^^f 
X2 = {l + X- ^/l + x'^f 

X3 = 1 
X4 = 1 



1 



m = ^^i,/3i,o,o) 

n2 = :^(l,/32,0,0) 
U3 = (0,0,1,0) 
M4 = (0,0,0,l) , 



(A.^ 
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where 



X 
/92 



7? 
7i 



1 + /3? 
1 + 



(A.9) 



A. 2 Weak ferromagnet 

In order to get a symmetric transfer matrix, we insert an identity matrix 
between the single state matrices m in (| 

1 ^ 




1 



2x2 



1 / \ 
This leads us to a modified single state matrix 

~ I |1) sgn a^ 

]a||T) 



(A.IO) 



m 




which obviously produces the same global state (|2] 
matrix for m is 





/ 


1 


a 













o 


1 




















0" 
























Its eigensystem 


is 










XI = 1 + 


a 






1 


(1,1,0,0) 


X2 = 1 




a 


? 


U2 




(1,-1,0,0) 


X3 = 0- 




; 


U3 


= (0,0,1,0) 


X4 = 0- 




; 


U4 


= (0,0,0,1) . 



(A.ll) 

as m. The transfer 
(A.12) 



(A.13) 



A. 3 Dimerized antiferromagnet 

As G and M alternate in the global state (|38|), one can use the transfer 
matrix for M • G or for G ■ M. However, since the transfer matrix for M ■ G 
would be a 9 X 9-matrix, it is more economical to use the transfer matrix 
for G ■ M, which is a symmetric 4 x 4-matrix, namely 



V 



2a2 + 4 



2a2 + 4 + 5 

















4/ 



(A.14) 



It has the following eigenvalues and eigenvectors: 
Xi = (a2 + l)2 + 8 , ^ii = ^(1,1,0,0) 



X2 

X3 
X4 



(a^ 
4 
4 



1 



U2 



1,0,0) 
(0,0,1,0) 
(0,0,0,1) . 



(A.15) 
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